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ABSTRACT 

We show that the high optical depth of the intergalactic medium to Lya photons be- 
fore the Epoch of Reionization results in a negligible magnitude for the Wouthuysen- 
Field effect produced by a radiation source on its distant surroundings, unless (1) 
the scattering medium has sufficient time for the impinging resonance line photons to 
establish a steady-state frequency distribution or (2) the scattering gas is undergo- 
ing internal expansion or has a peculiar motion of tens to hundreds of kms~^ away 
from the source. Because of intergalactic attenuation, discrete structures will receive 
only radiation from a source displaced from the resonance line frequency by typically 
hundreds to thousands of Doppler widths. The incident radiation must diffuse across 
the resonance line to produce a substantial scattering rate. We present steady-state 
solutions in the radiative diffusion approximation for the radiation field trapped in 
a clump of gas and show that this may result in an enhancement of the strength of 
the Wouthuysen-Field effect by as much as a factor of 10^ over the free-streaming 
(single-scattering) limit. Solutions to the time-dependent diffusion equation, however, 
show that the timescales required to establish a steady state will generally exceed the 
lifetime of the sources, resulting in a substantially reduced scattering rate. In the pres- 
ence of internal expansion, a steady state may be established as photons are redshifted 
across the resonance line and into the red wing, and significant enhancement in the 
scattering rate over the free-streaming limit may again be produced. Alternatively, 
a substantial scattering rate may arise in systems with a peculiar motion away from 
the source that redshifts the received radiation into the resonance line centre. As a 
consequence, at epochs z < 30, when collisional decoupling of the hyperfine structure 
of hydrogen from the Cosmic Microwave Background is small except in dense regions, 
and prior to the establishment of any large-scale diffuse radiation field of resonance 
line photons, the 21cm signature from the Intergalactic Medium produced by the 
Wouthuysen-Field effect will in general trace the peculiar velocity field of the gas in 
addition to its density structure. 

Key words: atomic processess - cosmology: theory - line: formation - radiative 
transfer - radio lines: general - scattering 



1 INTRODUCTION 

Following the Recombination Era at z ~ 1100, the baryons 
produced in the Big Bang were largely neutral. By z > 6, 
the spectra of high redshift Quasi-Stellar Objects (QSOs) 
show that the hydrogen had become ionized IjBecker fc et al.l 
[iool). The history of structure formation in the wide red- 
shift expanse between these epochs, when there were few 
if any sources of radiation, is largely unknown. In princi- 
ple, 21cm emission from intergalactic hydrogen could reveal 
the evolution of structure in the baryons during these cos- 
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mic "Dark Ages." Because the baryons are cold, they would 
closely trace the evolution of the dark matter, so that 21cm 
imaging could tra ce the growth of struc t ure following the Re - 
combination Era (iHogan fc Reedl 19791 : [Scott fc Reej|l99Gl '). 
Because of strong coupling between the hyperfine structure 
of the hydrogen and the Cosmic Microwave Background 
(CMB), however, a mechanism that decouples the spin tem- 
perature of the hydrogen hyperfine structure from the CMB 
temperature must be active, otherwise the hydrogen is ren- 
dered invisible: it absorbs and re-emits the CMB radiation 
at the same rate, leaving it indistinguishable from the CMB. 
At redshifts 2 > 30, collisions between hydrogen and elec- 
trons and other hydrogen atoms are adequate to begin de- 
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coup l ing the spin temper ature from the CMB l|Scott fc ReesI 
ll990l : lMadau et aLlll997l 'l. Except in dense regions, however, 
at later times coUisional decoupling is inadequate. 

Once the first sources of radiation begin to turn on, the 
scattering of Lya photons by the neutral hydrogen offers an 
alternative means of decoupling, through the Wouthu ysen- 
Field efi'ect (WFE) l|Wouthuvsenl Il952l : iFieldl 1 19581 ). As 
sources turn on in sufficient number to begin reionizing the 
Intergalactic Medium (IGM) , they will produce a combined 
intensity of Lya photons sufficient for coupling the spin 
temperature to the light temperature of the Lya photons 
rather than the CMB, rendering the intergalactic hydrogen 
visible against the CMB in either emission or absorption 
l|Madau et al.ll 19971 ). 

The discovery of the End of the Dark Ages (EDA) and 
the onset of the Epoch of Reionization (EoR) has become 
one of the paramount goals of a new generation of radio 
telescopes, such as the LOw Frequency Array (L0FARj3, 
the Murchison Widefield Array (MWA) 0, the Primeval 
Structure Telescope/21 Centimeter Array (PaST/21CMA) 
0, the Precision Array to Probe EoR (PAPERjl and 
a possible Square Kilometre Array (SK Afl Reviews of 
this rapidly growing are a are provided bylLoeb fc Barkanal 
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Pa that drives the WFE as the integrated intensity from 
cosmologically distributed sources, assuming that photons 
blueward of Lya emitted from a source will contribute 
their full amount to Pa once they redshift into the local 
Lya resonance. In fact, this provides a lower limit (as- 
suming photons are not destroyed by, eg, dust absorp- 
tion) . The multiple scattering of resonance line photons will 
in general enhan c e the radiation field. For pure Doppler 
scattering, iFieldl I 1959al ) argued the rate is enhanced by 
the number of scatters a photon undergoes before being 
randomly scattered sufficiently redward of line centre to 
escape. The estimate neglected scattering in the Lorentz 
wings, however. It also assumed that the IGM takes part 
in the homogeneous and isotropic expansion of the Uni- 
verse. In fact the IGM is clumpy, with structures break- 
ing away from the cosmological expansion. Efforts are un- 
derway to estimate the radiation field and scattering rate 
in a cosmological context through Monte Carlo simulations 
of the scattering of resonance lino photons in an inhomo- 
geneous medium (Chuzh oy & Zheng 20o3; IPierleoni et al.l 
l2007l : ISemelin et alJl2007D . In this paper, we use analytical 
means to explore some of the consequences of intergalactic 
structure formation for the WFE as a means of generating 
an intergalactic 21cm signature. 

In the next section, we discuss the optical depth of the 
IGM to resonance line photons and the implications for the 
scattering rate. Approximate steady-state solutions to the 
radiative transfer equation are derived in § [S] and time- 
dependent solutions in § m A discussion of the solutions 
and applications are provided in § O along with a summary 
of our conclusions. 
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Figure 1. (Upper panel) The optical depth of photons emitted 
by a source at redshift zg and received at the frequency u near the 
resonance line frequency of Lya at redshift z. The frequency 
is expressed a,s x = {u — vq)I i^vj^,. The curves are for T = 10 K 
(solid), 100 K (dashed) and 1000 K (dot-dashed). Also shown is 
the approximation ~ x\lx for the T = 10 K model (dotted 
line, nearly coinciding with the solid line). (Lower panel) The 
value of X at which Ti, = \ for zs — ^ < z < zs- The curves are 
labelled as in the upper panel. 



2 THE WFE IN THE INTERGALACTIC 
MEDIUM 

The optical depth through a homogeneous and isotropic 
expanding IGM of a photon emitted by a source at red- 
shift zs and received at redshift z at freque ncy v > vq , 
where vq is the resona nce line frequency, is l|Fieldl Il959al : 
iGunn fc Petersonll 19651 ) 
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I dip 
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where ni (z') is the number density of scattering atoms in the 
lower level at epoch z' ,a — ire^ fiu/{mec) ~ 0.0110 cm^ Hz is 
the total resonance line cross section, where /;„ ~ 0.4162 is 
the upward oscillator strength for hydrogen Lya, ipv{a,v) 
is the Voigt line profile normalized to J dvipv{cL,u) — 1, 
a ~ 0.0472r"^/^ is the ratio of the decay rate to the 
Doppler width Ai/d ~ vob/c, where h = {2k-BT /m-aY^^ 
is the Doppler parameter for hydrogen gas at temperature 
T and c is the speed of light, and Ip is the proper path 
length. In the Lorentz wing, expressed as a function of 
X = {v—vo)/Avd, the dimensionless Voigt profile (f)v(a, x) = 
(Ai/i5)((3v(a, ly) is well-approximated as (l)v{a,x) ~ a/ijix^). 
The differential proper line element evolves according to 
dlp/dz ~ {c/ Hq)Q,^^'^ (1 + z)~^^'^ in a flat universe at red- 
shifts for which S7m(l + z)^ dominates the contribution from 
the vacuum energy, where f2m is the ratio of the total mass 
density to the critical Einstein-deSitter density, and where 
Hi) — lOO/ikms"^ Mpc~^ is the Hubble constant today. 
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Photons emitted by the source at a frequency Ve > 
1/0(1 + zs) / (1 + z) will scatter in the blue Lorentz wing (ex- 
cept for those for which the inequality is nearly an equality, 
in which case they will scatter through the Doppler core). 
In the blue wing the optical depth due to Lya scattering 
by neutral hydrogen is 



1.75 X 10' 



cniiO) (l + zf/'' 



1/2 



y 



1/2 



,1/2 



+ 2 log 



(l-yi/2)(i + ^i/2) 



(2) 



where y = v/vo and u = y{l + zs) / (1 + z), and all quantities 
are assumed to be in cgs units. For < xb/c « {zs — z)/{l+ 
z) « 1, the expression simplifies to 
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for a mean line-of-sight IGM temperature Tigm for a uni- 
verse with Qrn = . 3 and baryon density Qth^ ~ 0.022 
l|0'Meara et al.ll2006l : [Spergel fc et al.ll2007i ). As an illustra- 
tion, Tu{x) for a source at = 8 is shown at z = 7 in Fig.[T] 
(upper panel). The optical depth is extremely high until u 
is well displaced from the line centre. The values of xi, at 
which T,j — 1, are shown in the lower panel of Fig. [T]for a 
range of redshifts. 

Photons emitted by the source with frequencies between 
vo < < i'o{l + Zs) / {1 + z) will pass through the resonance 
line frequency en route to the gas at z. Upon passing through 
the resonance line, the photons will essentially be completely 
scattered out of the line of sight. As a consequence, essen- 
tially no photons will be received in the frequency range 
!/o(l + + zs) < V < vo (the Gunn-Peterson effect). A 

general expression for the optical depth through the Doppler 
core is provided in Appendix 1X1 

Photons emitted sufficiently redward of uq to avoid the 
Doppler core may still be scattered in the Lorentz red wing 
of the surrounding hydrogen. This will produce a halo of 
scattered Lya photons around a central source like a QSO 
llLoeb &: Rvbickil 1 19991) . For f requencies a; < near line 
centre, iLoeb fc Rvbickil l| 19991 ) show that the photons will 
build up an energy density as they diffuse through the sur- 
rounding IGM of the form u^, ~ (_a;:)-9/2g-['-/'-c(a^)] ^ ^here 
r^{x) = (-a;)^''^(2/3)[(6/i?)An,fp(a;)]^/2 scales like the har- 
monic mean between the distance over which Hubble expan- 
sion produces a velocity difference matching the Doppler pa- 
rameter b and the scattering mean free path of a resonance 
line photon in the Lorentz wing 



tvAud 2 
Amfp(a;) ~ X , 
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where n; is the number density of atoms in the lower level. 
The energy density is exponentially suppressed for a; — > 0~ 
at a fixed position, resulting in a negligible Lya scattering 
rate except very near the QSO. Numerically, 



rc{x) ~ Q.6pc{-xf^'^h 
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so that a substantial Lya scattering rate is confined to a 
small region around the source. The region is so small in 
fact, that the assumption that the surrounding IGM takes 
part in the Hubble flow will not be valid. More likely the gas 
will have been ionized and violently disturbed by the QSO 
itself. 

Finally, photons emitted sufficiently redward of vq at 
the source to escape scattering in the Lorentz wing will pass 
freely through the IGM, where they will be received at red- 
shift z at frequencies x < —xo, where 



Xo 



+ zs 
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The vast majority of photons from the source emitted 
blueward of, but near, the resonance line frequency never 
arrive near the line centre, having been scattered out of the 
path along the way. The resonance line radiation imping- 
ing on the absorption sit^fl from the source scatters in the 
Lya resonance line at the rate per neutral atom in the lower 
state 



Pi = ac I dfipvia,!/)-^-^, 
Jo "■P'^ 



(8) 



where ftp is Planck's constant and is the local specific 
energy density of photons received at z from the source at 
Zs with specific luminosity 
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where tl is the luminosity distance between the source and 
the gas at redshift z. The scattering rate assuming none of 
the photons have been scattered out of the line of sight en 
route to the absorption site (the free-streaming limit) is 



Pi{t^ =0) -o- 



Lv/hpv 
^itr\ 
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The ratio of the scattering rate for ~ xi/x to the free- 
streaming value is then given by 

P;(t^) ^ dxa-K~''^x~^ exp{-xi/x) 

Pi{ti, = 0) S^oo dx(j>v{a,x) 

~ 7.0xl0-^ft(l + ^)"'^'(^)''', (11) 

where Ta is the temperature of the absorbing gas. The low 
value poses a fundamental problem to the effectiveness of 
the WFE as a means of decoupling the spin state of the gas 
from statistical equilibrium with the CMB. Only if photons 
have sufficient time to diffuse across the line centre will the 
scattering of Lya photons be an effective means of decoupling 
the spin temperature from the CMB in comoving objects. 



3 THE WFE IN DISCRETE OBJECTS 

3.1 Steady-state solutions without atomic recoil 

The assumption of a homogeneous and isotropic uni- 
verse, while a fair approximation over large scales, breaks 

^ Strictly speaking, the photons are scattered. The region at 
which the scattering rate is computed is referred to for simplicity 
as the absorption site or the absorber. 
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down on small scales: the IGM is clump y at the redshifts 
from which the 21cm s ignal is expected l|Tozzi et al.ll200ol : 
iGnedin fc Shaved l2004l l. For a sufficiently short mean free 
path, photons will be trapped within a discrete system and 
diffuse in frequency and space as they scatter within the 
system. In this section, we consider the energy distribution 
the radiation will reach given adequate time to relax to a 
steady state. 

For a Lya photon to become trapped within an ab- 
sorber of lengthscale La, the mean free path must mini- 
mally satisfy Amfp < La, for which r,^ = La /Xmip > 1- 
As the optical depth in the wing may be expressed as 
r,j — Tv~^^^{aro)x~^, where tq = niLaTv'^^'^(j{Auo)~^ is the 
optical depth at the line centre, this c orresponds to pho - 
tons escaping with x > 7r~^''*(aro)^^^ (jOsterbrockl 1 1962 1 . 
Radiative transfer solutions for sources embedded within a 
slab suggest this may be somewhat too restrictive for sys- 
tems with very high line centre optical depths, in which case 
the photons escape through spatial diffusion rather than 
frequency diffusion. Steady-state solutions in the diffusion 
approximation for scattering in the Lorentz wings in very 
high optical depth systems show that photons escape the 
slabs with frequencies typically on the order of (aro)^^'^ 
l|Adams|[l973 : iHarringtonl Il973| - iNeufeldl Il990l ). An escape 
frequency of Xesc = a;«(aTo)^''^ corresponds to a mean free 
path smaller than the thickness of the slab by the factor 



/ 



.1/2 



-1/3 



OmSxlT^^N'^^^, where Ta is the 



temperature of the absorber and Nig is the column den- 
sity in units of 10^^ cm~^. (This corresponds to a region 
24 kpc across for the mean intergalactic hydrogen density 
at 2; = 8, and smaller for overdense systems.) As / depends 
only weakly on Ta and A^hi, while j:* depends on aro and 
the radiation source parameters of any particular problem, 
we express the escape frequency simply as 

x^sc =i 395(10/)'/'AfiV'T'-''", (12) 

noting that / will in general vary with the depth within the 
absorber, and must be solved for according to each particu- 
lar configuration, but will be of order 0.01 — 1 for applica- 
tions we consider. Comparison of Eq. (I12p with Eq. (|4]) for 
xi shows that photons with r,^ < 1 in the IGM will gener- 
ally not be trapped within a structure comoving with the 
expansion of the Universe unless the structure has a neutral 
hydrogen column density A^ig > A''i9,crit, where 



A-i9,„it ==; 0.29ft~'(10/)~'(l + zf (-^] 

\ 1 IGM / 



(13) 



If x\ < Xosc, the absorber will behave as a photon 
bucket, allowing the energy density of the radiation field 
at frequencies x < x^sc to build up with time as radiation 
from the source becomes trapped. As the photons scatter 
within the absorber, they will diffuse in frequency and space. 
Because photons with \x\ > Xcsc will escape, the energy den- 
sity of the radiation trapped within the absorber will reach 
a steady state when the rate of incoming photons balances 
the rate of escape. In principle, the energy density that de- 
velops could produce a substantial scattering rate, even ex- 
ceeding the rate estimated assuming no photon losses from 
the source due to scattering by the intervening IGM. 

The depth to which photons at frequency xi penetrate 
an absorber depends on the optical depth of the absorber. 
If it is optically thin at xi, the photons will stream through. 



If optically thick, the photons will scatter near the sur- 
face of the absorber. Photons blueward of xi, however, will 
penetrate more deeply, as the cross section diminishes like 
1/x'^. In general, photons will be injected at varying layers 
within the absorber, with bluer photons injected at increas- 
ing depths. 

We would like to estimate the evolution of the radia- 
tion field within an absorber. No solution to this problem 
exists in the literature. The full radiative transfer problem 
through the slab is an involved one, as the radiation field will 
vary with both the optical depth through the slab and with 
frequency. It is particularly important in our application to 
obtain the solution across the Doppler core to estimate the 
Lya scattering rate. Instead of solving the full problem, we 
seek approximate solutions in which we treat the scattered 
radiation field as locally isotropic, neglecting the spatial dif- 
fusion of the radiation through the absorber. We show in 
Appendix|B]that the resulting radiative transfer equation in 
the diff usion approximatio n is formally identical to the equa- 
tion of IHarringtonl (|l973l ) for the problem of a slab with a 
uinform distribution of spectrally flat sources. The solution 
we obtain neglecting the spatial diffusion of the radiation 
agrees well with Harrington's solution which includes both 
spatial and frequency diffusion. A detailed quantitative de- 
scription will require a more exact solution to the problem, 
but our approximate approach allows us to explore qualita- 
tively several effects of interest. We would like to estimate 
the energy density that may be achieved as a function of 
lesc- We also wish to examine the effect of internal motions 
on the radiation energy density across the line centre and 
the possible role of atomic recoil. Lastly, we would like to 
estimate the time it takes the radiation fleld to establish a 
steady state. All of these effects extend well beyond those 
explored by existing slab solutions. 

We estimate the radiation density that builds up within 
an absorber by assuming a steady state between the rate 
of photons injected into the absorber and the rate of dif- 
fusion of the photons across the escape frequency, impos- 
ing the boundary condition u{x) — for > a;osc, where 
u{x) — AvdUv We assume isotropic scattering within a ho- 
mogeneous absorber, and use the diffusion approximation 
to describe the resonance line scattering of the photons. In 
equilibrium. 



ld_ 

2dx 



D(a, X 



dn{x) 



dx 



-Six) 



(14) 



where n{x) — Auon^, = Uy/h-pu, S{x) = 

{AuB)^S{u)/{nica), where Slv) is a source function, and 
D{a, x) is the diffusion coefficient. The second order moment 
of the frequency redistribution function gives for the diffu- 
sion coefficien t D(a,x) = (j)v (a,x) + {l/3)d^ (j>v{a,x)/dx^ 
(|Rvb icki & dc irAntoniclll994l '). The derivation of the diffu- 
sion equation, however, does not conserve photon number to 
the order of the approximation. As a consequence, the dif- 
fusion coefficient is left ambiguous. iRvbicki fc deirAntoniol 
(jl994!) advocate adopting D{a,x) — (f>v{a,x) for its simplic- 
ity, and we shall do so here. 

The source S{v) represents the rate at which photons 
received from the source are incident on the absorber at 
X < Xcsc, where they are trapped and diffuse in frequency. 
It is given by the scattering emissivity 
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Figure 2. The steady-state distribution of Lya resonance line 
photons for a source function S{x) = Soi^ — Xinj) with Zinj = 100, 
allowing photons with \x\ > Xqbc = 200 to escape, as given by the 
diffusion approximation, Eq. I I14I I. with D(a,x) = 0\/(a, x). An 
absorber temperature = 100 K is assumed. 
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dv' R{v' , i')n^i , 
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where R{i'' , v) is the photon redistri bution function for reso- 
nance line scattering dMihalaslllQTil ). Since scattering in the 
wings is coherent, R{v' , v) — tpv{a, i/')Sd{i'' — v), where 5d 
is the Dirac 5-function. Then 



S{v) 



nica<pv{ci, v)ni, 



1 L^/hvv 



X exp(— a;i/a;). (16) 



The dependence on x resuhs in a sharp peak at j: = xi/2, 
so that the source S{x) is well approximated by 



S{x) 



TIC xi ATir\ 

A5d{x - Xinj) 



Sd{x - 



.j) (17) 

where the normalization A — 2.5 x lQ^^h{l + 
z)~^/^T^J^[(^L„ / h-pv) / 'i-Kr\] (with all quantities expressed 
in cgs units), matches the source rate of Eq. p6|l integrated 
over all frequencies. Note that the dependence on the ab- 
sorber temperature cancels. Here, the frequency at which 
photons are injected into the absorber is represented gener- 
ally as Xi-nj to allow for any peculiar motion of the absorber, 
which could shift the peak of the received radiation either 
redward or blueward, and to allow for scattering within the 
absorber which will shift the peak of the radiation field pen- 
etrating to deeper layers further toward the blue. For a co- 
moving absorber, a;inj = x\/2 near the surface. The solution 
to Eq. (|14p using the wing approximation for <j)v{a,x) is 



n{x) 



3a ' 



^A 1 + 



< X < Xin 



•^ini X Xq 



The solution is shown in Fig. (2] We find that adopting 
D{a,x) = (j>v{a, x) + (l/3)ci^(?!)v(a, x)/dx^ produces a nearly 
identical solution, agreeing to within a fraction of a per- 
cent for Xinj well out of the core. Numerically integrating 
the diffusion equation using the full Voigt profile produces 
a solution that also differs negligibly from Eq. (|18p . 

The ratio of the collision rate Pi, given by using Eq. (|18|) 



Figure 3. The steady-state distribution of Lya resonance line 
photons for a source function S{x) = S o {x — xi^ij) with a'inj = 100 
in an absorber with temperature Ta = 100 K, allowing photons 
with I a; I > Xesc = 200 to escape and including the effects of inter- 
nal motions with a dimensionless velocity gradient 7 = 1.7 X 10~^ 
(solid line), representing expansion. Virtually all the injected pho- 
tons are redshifted, escaping for x < —Xcbc- The solution for 
Xinj = 25 and x^sc = 50 for 7 = 1.7 X lO"'^ is also shown (long- 
dashed line) , showing a weak tail at x > Xi^j. The corresponding 
case with contraction, 7 = — 1.7x lO"'^ (short-dashed line) shows 
a more complex behaviour as contraction curtails the redward 
diffusion. 

in Eq. ([8]), to the free-streaming collision rate Pi{0) of 
Eq. dlO} is then 
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1 a 

Xl TvA 
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-1 3/2 



(1 + z)Ta 



1/2 



*'CSC 



(|AT|| below for 
assumed, and 



where we have used the approximation Eq. 
the Voigt profile with a;csc > a;inj 3> Xm 
neglected the terms following the leading. The enhance- 
ment may be substantial. For A'^ig = A^ig.crit, Pi,ss/Pi{0) — 
1.5 X 10*/i"^(l -f zy-^T{^^^{l - xf^j/x'i^^), so that even for 
a warm IGM temperature of Tigm = 1000 K a substantial 
boost in the scattering rate will result. 

The solution Eq. (|18|) assumes velocities internal to the 
cloud are negligible. To estimate the effect internal motion 
may have on the frequency distribution of the photons, we 
approximate the motion as an isotropic expansion (or con- 
traction) within the cloud characterised by a uniform veloc- 
ity gradient Ha (which in general will not be the same as the 
(18) Hubble param eter). In the diffusion approxim ation, Eq. (|14|l 
is modified to l|Rvbicki fc deU'Antoniolll994l ') 
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2d^ 



2'yn{x) + D{a 



dn{a 



dx 



= -S{x), 



(20) 



where 7 = uoHa/(unic) is (c/6) times the ratio of the scat- 
tering time of the photons at line centre to the expansion 
time of the gas. For this solution, the steady state 

applies to the comoving number densi ty of photons per 
frequency l|Rvbicki fc deirAntonici|ll994 ). Overdense struc- 
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Figure 4. The steady-state distribution of Lyo resonance line 
photons in the diffusion approximation for a source function 
S{x) = SdIx — Xinj) with Xinj = 100, allowing photons with 
\x\ > XcBc = 200 to escape and including the effect of atomic 
recoils. Shown are the solution with e = 0.0079, corresponding 
to Ta = 10 K (solid line), the solution with e = (long-dashed 
line), and the first order solution no{x) + eni{x) for e = 0.0079 
(short-dashed line). 



tures are typically contracting, but the internal motions 
along a filament m ay expand along the filament in regions 
l|Zhang et al.lll99^ . For an overdensity ten times the cosmic 
mean at z = 8, an internal velocity gradient of lOOkms"^ 
over 100 kpc corresponds to 7 ~ 1.7 x 10"^. Contraction by 
the same magnitude corresponds to 7 ~ —1.7 x 10~^. 
The solution for the source ^(a::) = A5d{x — Xinj) is 



n{x) 



9-1 {x, 



^CSC 



) 5^7 (-^i "^inj ) ) 3^inj X <^ X 

)f'r{x,y)/ fiix. 



CSC , 



,(21) 
,y) and 



'—-{x'^ — y^)]. The solution is extremely 



where g-t{x,y) = {A/y)fj(x, 
fj{x,y) ^ l~exp^ 3^ 

sensitive to "fx^^^/a. For 7 ^ 7crit = Sa/{2-KX^^^), a broad 
wing of amplitude A/"f redward of Xinj will result as photons 
are redshifted across the line centre. The distribution will cut 
off sharply at 2; > Xinj. For an absorber with Ta = 100 K 
(a ~ 0.00472), and a;csc = 200, 7crit ^ 2.8 x 10"^°. The so- 
lutions for 7 = 1.7 X 10"^ with (a;inj, a;csc) = (100,200) and 
(25, 50) are shown in Fig. O A substantial decrease in n{x) 
is produced compared with the 7 = case. A significant en- 
hancement of the Lya collision rate over the free-streaming 
limit, however, may still result, with Pi^bb/Pi{0) — a/^n'yx-i). 
For 7 < 0, a broad wing blueward of Xinj will form as pho- 
tons are blueshifted in the contracting gas. A complex pro- 
file, however, may form redward of redward diffusion 
is resisted by blueshifting. Such a case is illustrated in Fig. [3] 
for 7 = —1.7 X 10~^ with a;inj = 25 and a:esc = 50. 



3.2 Steady-state solution with atomic recoil 

The above solutions neglect the effect of atomic recoil on 
the evolution of the radiation field. Since the photon mo- 
mentum (~ hpuo/c) is small compared with the momentum 
of the hydrogen atoms (~ mub), the effect should be small, 
but may result in a significant heating rate when the gas is 
stiU cold dMadau et al.lll997l : IChen fc Miralda-Escudj|2004 
lMeiksinll2006l ). Allowing for atomic recoils modifies Eq. I|14p 
to 



ld_ 

2dx 



2e(f>v{a, x)n{x) + D{a, x 



dn{x) 



dx 



= -S{x), 



(22) 



where the recoil parameter e = hpU()/(2kBTmiic^)^^^ ~ 
0.025r~^''^ is the ratio of the momentum of a resonance line 
photon to th e momentum of an atom mov i ng at the ther - 
mal velocity l|Rvbicki fc deirAntoniolll994 lMeiksinll2006l ). 
Eq. (|22|l expresses only the first order correction in e. The 
solution for a source S{x) = A5d{x — Xinj) is 



i{x) = 



X\y].\ X, 



CSC ) 



(23) 



where ge{x,y) = A[tv /{2e^a)]exp{~2ex)[f^{x) - 
/E(y)][/E(a;osc) - /£(a;inj)]/[/,(a;csc) - .My)], and 
f^{x) = (1 - 2ex + 2e'^x'^) exp{2ex). Here, the wing 
approximation for (pv{a,x) is adopted, but we find a nearly 
identical result allowing for the Doppler core as well. 

The solution for xinj = 100, Xcsc = 200 and e = 0.0079, 
corresponding to T ~ 10 K, is shown in Fig. The re- 
sult differs substantially from the solution for e = 0. In 
fact the solution cannot be trusted, as it is not consis- 
tent with the original order of Eq. (|22|l in e. This may 
be demonstrated by adding to 2e<l)v{a, x)n{x) the quan- 
tity e^y{x), where y{x) is an unspecified function represent- 
ing the higher order corrections to the recoil, and which 
in general may also depend on n{x). Inserting the series 
solution n{x) = "Sj^o ni(x) into Eq. (|22p and equating 
equal orders in e recovers naix) as given by Eq. (|18|l . and 
produces the first order equation {d/dx)[2())v{a,x)no{x) + 
D{a, x)dni{x)/dx] — and the second order equation 
{d/dx)[2(j)v{a, x)ni{x) + y{x) + D{a, x)dn2{x) /dx] — 0. The 
first order equation gives ni{x) in terms of no{x), subject 
to the boundary conditions ni(— Xesc) = ni{xasc) = 0. The 
second order equation shows that n2{x) depends on the un- 
specified (neglected) function y{x). Thus the solution is valid 
only up to n{x) = no (a;) -I- eni{x). If higher order terms are 
substantial, then the solution for n{x) cannot be trusted. 
The solution no{x) + eni{x) for the example above differs 
substantially from the solution n{x), and indeed is not even 
positive-definite, as shown in Fig. 2] demonstrating the so- 
lution is not valid. 

Including higher order terms involves a non-trivial ex- 
pression for the redistribution function including non-linea r 
terms in the recoil parameter (lBaskdri98ll : [Meiksinll2006l ). 
(Additional higher order terms in the recoil parameter arise 
from relativistic corrections, but these are reduced by factors 
of b/c.) For temperatures T > 100 K, the discrepancy be- 
tween the full solution and the first order solution is smaller 
than about 20% for Xinj = 100 and x^sc = 200. In general, 
the correction will be small for £2;osc ^ 1, and we caution 
that if this is violated the steady-state solutions may not be 
valid. We do not pursue further consequences of the atomic 
recoil solution here. 



3.3 The role of bulk peculiar velocities 

Absorption systems with a bulk peculiar velocity compo- 
nent away from the source will see the photons redshifted 
into line centre. In this case, the scattering rate of resonance 
line photons may be appreciable even in an optically thin 
absorber {Nig < A^ig^crit), so that the build-up of a strong 
radiation field at line centre through photon capture and 
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Figure 5. Surface of constant Lyo scattering rate Pi produced by 
a continuum source to the left witliin a collapsing halo of neutral 
hydrogen. The inner surface terminates at the accretion shock. 
The circle represents the turnaround radius rta of the halo, with 
the solid line representing the hemisphere of gas with peculiar 
velocity redshifted relative to the source and the dashed line the 
hemisphere with blueshifted peculiar velocity. The surface shown 
corresponds to a maximum radius of 0.59rta. The figure is tilted 
by 5° to reveal the 3D structure of the surface. 



diffusion is not necessary for significant WFE decoupling of 
the spin state from the CMB. The peculiar velocity required 
to bring the absorber into the peak of the source function 
S{i') in Eq. (|16|l is given by t;poc — ^xiba, where x\ is given 
by Eq. (|4]). For pure thermal broadening (assuming negli- 
gible internal velocity structure within the absorber), this 
corresponds to 



1/2 



Uh-\l + zf/^ ( ^\ kms- 

V TlGM / 



(24) 



This is comparable to the typical peculiar velocities of non- 
linear cosmological structures at the epochs of interest. As 
a result, the WFE will produce a 21cm signature that maps 
the peculiar velocity field of the structures within the IGM. 

As an illustration, we consider the signal produced by a 
collapsing halo, such as around a forming galaxy or galaxy 
cluster, illuminated by an external radiation source, as de- 
picted in Fig. [5] 

The collision rate from a source of specific luminosity 
Li, on an absorber with Doppler parameter ba at luminos- 
ity distance moving away from the source with peculiar 
velocity v is 



hp V A-Krj^ 



hpv 47rr? 



exp 



di/ifv \a,i'(l -)j exp(— r,y) 



vo cos{6) 



(25) 



using Eq. (|3]), approximating <j)v{a,x) ~ tt"^''^ exp(— 
noting that the resulting exponential in the integrand peaks 
sharply near x — v/ha assuming vo/ba — xi or larger, Taylor 
expanding the argument of the exponential to second order 
in a; — v/ba, and expressing the peculiar velocity as v = 



vqcos{0), the projected accretion velocity along the line of 
sight to the source. 

In general, «o will vary with distance from the centre of 
the collapsing halo. As a definite case, we approximate the 
peculiar ve locity of the accreti ng gas using the self-similar 
solution of iBertschingeij (|l985l ) for adiabatic accretion of a 
7 = 5/3 gas of negligible density onto a dark matter halo 
with a comoving centre in an Einstein-de Sitter cosmology. 
The velocity field in the region between the accretion shock 
and the turn-around radius at cosmological time t is well- 
approximated by 



i'o(r) 



t log^oA, 



■ log 10 ^ 



(26) 



where rta is the turn-around radius, A = r/rta. As = 
rs/rta — 0.3472, where is the radius of the accretion 
shock, and Vmin — —1.433 characterises the infiow velocity 
of the gas just before passing through the accretion shock. 

Contours of constant values of Pt correspond to values of 
constant projected accretion velocity. If at a radius r towards 
the source (6 — 0), the accretion velocity is a fraction / of 
its maximum infall velocity, then contours of fixed Pi will 
correspond to the surfaces 



A = A{ 



/soc(e) 



(27) 



extending over the angular range < ^ < arccos(/). 
The surface corresponding to f — 0.5 is shown in Fig. [5] 
with a maximum angle 6max = 60° and maximum radius 
Ao = A{ = 0.59. The loci of constant Pi correspond to arcing 
sheets. The exponential sensitivity of Pi to v for vo/ba — xi 
(see Eq. [55]) ensures a surface of narrow width will dom- 
inate the 21cm signal produced. Such a surface could be 
mistaken for the region of neutral hydrogen just outside an 
H 11 region surrouding a source of ionizing radiation. A vary- 
ing gas temperature, and therefore varying ba, in the accret- 
ing region will further complicate the signal. If vo/ba >> xi, 
then the entire solid hemisphere will be illuminated by the 
source (except for a thin circular wafer at 9 < 90°). 

A similar effect may apply to an expanding void. Since 
voids expand faster than the Hubble expansion, the peculiar 
velocity on the far side of a void from a distant source could 
redshift the received radiation into line centre. A shell of 
direct scattering may then result even without including the 
reddening effects of radiative transfer within the void. 



4 THE EVOLUTION OF THE WFE IN 
DISCRETE OBJECTS 

It was assumed for the solutions in § [3] that the photons 
are able to achieve a steady state as they scatter within the 
absorbing system. In this section, we compute the timescales 
required to establish a steady state by solving the time- 
dependent diffusion equation 



dn{x,T) 
dr 



2 dx 



D{a, x) 



dn{x, t) 



dx 



+ S{x) 



(28) 



where r = t/ts, and ts = Aun / (riiac) ~ 3.2nY^Ta^^ s is the 
scattering time of the resonance line photons at line centre. 

We solve Eq. (|28|) using a Crank-Nicholson scheme. We 
tested the scheme against the time-dependent solutions for a 
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Figure 6. (a) Convergence of the time-dependent solution to 
the diffusion equation for a source S{x) = 5]j{x — a'inj) with 
Xinj = 100, Xesc = 200 and Ta = 100 K. Shown are the steady- 
state photon frequency distribution, Eq. 1181 1 (solid line), and the 
photon frequency distributions at the times their values at z = 
take on 50% (dashed line), 90% (dotted line) and 97% (dot-dashed 
line) of the steady-state value, (b) The evolution of the ratio of 
the photon frequency distribution at line centre {x = 0) to the 
steady-state value, as a function of time in units of the scattering 
time at line centre. 




Figure 7. The convergence time in units of the scattering time 
at line centre for the photon frequency distribution with source 
S{x) = So{x — Xjnj) having x^^j = x^sc/'i to reach 50% (dashed 
line) and 90% (solid line) of the steady-state value at x = 0. A 
temperature of Ta = 100 K is assumed. 



constant 5- function source at line centre in an infinite homo- 
geneous medium and a flat continuum source i n the diffusion 
approximation ("Rvbick i fc dell' Antonio! 1 1994h . The results 
applied to the source S{x) = Snix—Xinj) with Xinj = 100 and 
boundary condition n{x) = for \x\ > a;csc = 200 are shown 
in the upper panel of Fig. |6l The photon frequency distribu- 
tion takes on the approximate shape of the final steady-state 
solution Eq. (|18p as it converges towards it. The evolution 
of the photon frequency distribution at line centre (x = 0) 
is shown in the lower panel of Fig. [G] For sources that have 
been active only a small fraction of the time required to 
achieve full convergence, the scattering rate may be negligi- 
bly small. 

Dimensional analysis of Eq. ([28} suggests the time to 
establish a steady state for photons injected at a; = a;inj 
is Tss ^ 2a;f„j/D(a, Xinj) ~ (27r/a)a;f„j = {■K/8a)xt^^ ~ 
8.3Ta^^xt^^ for a^inj — a^esc/2. The time for the photon dis- 
tribution to converge to the steady-state value at a:: = 
found from the numerical integration of Eq. (|28|l is shown 
in Fig. [7| for convergence to 50% and 90% of the steady- 
state value, for a source S{x) — Sd{x — Xinj) with Xinj = 
Xf.sc/2 and Ta = 100 K. The convergence times scale like 
Tss,5Q ^ 210xt, (50%) and r,,,9o ~ 750xi, (90%), i 
agreement with the expected behaviour. 



5 DISCUSSION AND CONCLUSIONS 

The role played by the Wouthuysen-Field effect in decou- 
pling the spin temperature of the neutral hydrogen from the 
CMB temperature depends most crucially on the scattering 
rate of Lya photons. The spin temperature Ts is gener- 
ally a weighted mean of the colour temperature Ta, kinetic 
temperature Tk of the gas, and brightness temperature Tr 
of any incident radiation (such as the CMB) at the 21cm 
frequency, 



Ts = 



Tr -f ygTc, + VcTk 
1 + + 2/c 



where 



Pio n 

AlO Ta 



and 



Cio T* 
A^T^ 



(29) 



(30) 



are the Lyg and c ollisional pumping efficiencies, respec- 
3ly dF 

'^^ is the spontaneous decay rate of the 21cm tran- 



2.85 X 



umpr 

lively (|Fieldl 1 19581 : iMadau et all 1 19971 ). Aio 

10-15 s- 

sition, and T, = hpVio/kB, where vio is the frequency 
of the 21cm transition. Here Pio = 4P;/27 is the indi- 
rect de-excitation rate of the triplet hyperfine state in- 
duced by Lya photon scattering. The colour temperature 
of the radiation field is the harmonic mean temperature 
Ta = l/(T-i(i/)), where T„{u) = -~{hp/kB){dlnu^/dv)-\ 
weighted by Ui,Lpv(a , v) for an energy density of reso- 
nance line photons (jMeiksinl t2.006) . The coefficient Cio is 
the de-excitation rate of the triplet hyperfine state induced 
by atomic collisions. 

For a Lya scattering rate Pi exceeding the critical ther- 
malization rate 



27^ioTct 



: 6.8 X 10-1^ s" 



" V 9 " ^''^ 

where Tcmb ^ 2.73(1 + z) (|Mather fc et al.lll994l ) is the 
CMB temperature, the spin temperature will be driven to 
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the colour te mperature in the ab sence of strong coUisional 
de-excitation l|Madau et al.|[l997l ). The radiation field will 
typically rapidly thermalize with the neutral hydrogen, re- 
sulting in Tc Tk (|Field .f959b, : ,Meiksin 20061. The neu- 
tral hydrogen will then appear in either absorption or emis- 
sion against the CMB, depending on whether Tk is below 
or above Tcmb, respectively. 

Before the Epoch of Reionization, when the first radi- 
ation sources turn on at the End of the Dark Ages, regions 
around the sources may become visible through their 21cm 
signature. Because the photons produced by a source, such 
as a QSO or star-forming galaxy, near the Lya resonance 
will be scattered out of the line of sight over short distances, 
the scattering rate of Lya photons will be much smaller than 
Pth except extremely near the source. Only photons emit- 
ted sufficiently blueward of the Lya resonance frequency 
will survive to any distant structures, where they will be 
received well in the blue Lorentz wing. If a structure has a 
sufficient optical depth in the wing to the photons received, 
it will trap them and rescatterings within the structure will 
allow the photons to diffuse in frequency across the reso- 
nance line. The radiation field may be built up in this way 
until it reaches a steady state in which the rate of incoming 
photons is balanced by the rate at which photons diffuse suf- 
ficiently far from the resonance line frequency to escape in 
the wings where the structure is optically thin. The resulting 
photon scattering rate in the diffusion approximation may 
greatly exceed the estimate Pi{0) in the free-streaming limit, 
for which it is assumed no resonance line photons emitted 
by the source have been scattered out of the line of sight 
before reaching the absorber. 

We identify several factors, however, that may reduce 
the energy density of the resonance line photons from 
the steady-state diffusion approximation estimate for static 
structures: 

(1) Internal expansion within the structure, charac- 
terised by the dimensionless expansion parameter 7, will 
redshift the photons into a red wing, where they will es- 
cape, on a timescale much faster than the diffusion time, 
resulting in a much reduced energy density. In the diffusion 
approximation, the scattering rate will be reduced compared 
with the 7 = case by the factor (3a/7r7)/(a;|?gt. — a;f„j) for 
7 ^ 7crit = 3a/(27r2;esc), where Xesc is the dimensionless es- 
cape frequency at which photons escape the structure, and 
Xinj is the frequency at which photons are predominantly 
received from the source. Contraction within the structure 
may altogether prevent photons from diffusing across the 
resonance line frequency. 

(2) The scattering rate may be reduced by atomic recoil. 
Estimating the magnitude of this effect, however, generally 
requires including non-trivial higher order terms in the recoil 
parameter e than first order in the steady-state equation. 
We caution against introducing atomic recoils into Monte 
Carlo calculations and other steady-state solutions, as the 
results obtained may not be physically self-consistent if the 
solutions differ substantially from those obtained for e = 0. 

(3) The time scale to establish a steady state will typ- 
ically exceed the duration of high star formation rates in 
galaxies or the lifetime of QSOs. As a consequence, the col- 
lision rate will be reduced by factors of tens to thousands 
compared with the steady-state value. 



As an illustration, consider the region around a bright 
QSO with specific luminosity at the Lyman edge of 
10'^^ ergs s~^ Hz^^ pre-heated by soft x-rays ahead of the 
ionization front to a mea n temperature of T igm — 1800 K 
at a distance of 10 Mpc (|Madau et al.|[l997h . According to 
Eq. (gl, ~ 200 at 2 = 8 (assuming h = 0.7). A struc- 
ture 250 kpc across with a moderate overdensity of 3 and 
internal temperature of Ta ~ 120 K will have a neutral hy- 
drogen column density of A^hi — 3 x 10^" cm~^ and escape 
frequency Xcsc — 200, from Eq. (|12|l . For a power law spec- 
trum f"^'^, Eq. (|10|l gives for the Lya collision rate in the 
free-streaming limit Pi{0) ~ 8.7 x 10~^^ s~^, about an order 
of magnitude smaller than the thermalization rate Pth- The 
steady-state value from Eq. ((19} then gives Pi^bb — 10~* s~^, 
well in excess of Pth. The time to establish a steady state, 
however, is long. The scattering time of Lya photons in the 
structure is about ts — 8.5 x 10* s. According to Fig. [T] the 
time for the radiation field to achieve a steady state is then 
about 10^^ s, or 3 Gyr, greatly exceeding the expected life- 
time of the QSO. If the QSO had illuminated the structure 
for as long as 0.1 Gyr, according to Fig. [S] the steady-state 
value must be reduced by a factor of 0.05. The resulting scat- 
tering rate will still exceed Pth. If the system has a peculiar 
velocity away from the QSO, the rate could be quite differ- 
ent. Eq. ((24} shows that a peculiar velocity of 140kms~^, a 
plausible value, would shift the peak of the radiation from 
the source received by the structure to the resonance line 
frequency. The radiation field would then rapidly reach its 
steady-state value, and Pi^ss would greatly exceed Pth. 

A much more overdense structure at the same distance 
from the QSO would be able to achieve the steady-state 
scattering rate. For example, a virialized structure 35 kpc 
across with an overdensity of 200 and temperature Ta — 
1000 K would have a neutral hydrogen column density of 
3 X lO'^^ cm~'^ and a;csc — 220. The time for the radiation 
field to reach a steady state would then be about 0.2 Gyr, a 
plausible lifetime for the QSO. 

By contrast, it may be difficult to achieve a steady state 
in a void. For an IGM temperature Tigm — 100 K, xi = 825 
at z = 8. A region 3 Mpc across with a temperature of 
Ta = 10 K, and underdense by a factor of 3 will have a neu- 
tral hydrogen column density of 4 x lO'^" cm~^, a;csc — 800 
and a scattering time at line centre of ts — 2.2 x 10"" s. The 
time to reach a steady state is then 2100 Gyr. Even if the 
void is illuminated for as long as the age of the Universe, 
the radiation field will be negligibly small compared with 
its steady-state value. Since the void will expand, however, 
the i ncident radiation will be redshifted across the resonance 
line. iRvbicki fc deU'Antonid (|l99i) estimate a characteris- 
tic time for redshifting to dominate diffusion in an infinite 
homogeneous medium to be r-y ~ {af-y'^Y^^. Voids expand 
somewhat faster than the mean Hubble expansion. Adopting 
the Hubble constant at 2 = 8 for /i = 0.7, 7 ~ 5 x lO"*^. A 
steady state should be established for t ^ Tjts — 2 x 10"* yr, 
and the steady-state scattering rate will then match the free- 
streaming value, with Pi,ss/P(0) ~ 0/(717x1) ~ 1, which 
may be adequate for producing a 21cm signature in the pres- 
ence of sufficient sources. 

The examples above show that the effect of peculiar 
motions may lead to an enhancement in the Lya scattering 
rate over the free-streaming limit and produce the required 
decoupling of Ts from Tomb. Even if the scattering rate is 
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inadequate for inducing decoupling throughout most of the 
structure, it may be adequate along narrow loci within which 
the velocity field either permits the resonance line photons 
to accumulate near the resonance line frequency or shifts the 
resonance line frequency of the absorber to matching the in- 
cident radiation field from the source as filtered through the 
IGM. We presented an example of an accreting halo which 
produces a 21cm signature along a curved surface for which 
the impinging photons are received at the resonance line fre- 
quency due to the infall velocity of the accreting gas. Such a 
curved structure could be mistaken for a rim of neutral gas 
outside an H ii region in radio maps of the IGM. Similar 
curved surfaces may form on the far side of voids from a 
distant source. 

As a consequence of the above efi'ects, the 21cm signa- 
ture of the IGM will not simply trace the density structure 
of the IGM, but will trace the peculiar velocity structure 
as well. We have not considered multiple sources. In the 
presence of a large number of nearby sources, a greater frac- 
tion of the gas will be able to match the incident radiation 
field through their peculiar motions. In addition, the escap- 
ing radiation from the absorbers will contribute to the am- 
bient radiation field. In due course, a metagalactic diffuse 
radiation field may grow sufficiently strong to ensure the 
WFE is capable of decoupling the spin temperature from 
the CMB temperature throughout the IGM, and so ensure 
the production of a 21cm signature everywhere. The effects 
described here suggest that in the early stages, before any 
such metagalactic field has been established (and it has yet 
to be demonstrated that such a field will be established), the 
21cm signatures produced at the End of the Dark Ages, and 
possibly extending into the Epoch of Reionization, may only 
be interpreted with a thorough understanding of the evolu- 
tion of the energy density of resonance line photons within 
a clumpy medium, including the effects of peculiar motions. 



APPENDIX A: OPTICAL DEPTH THROUGH 
DOPPLER CORE 

Radiation from a source at zs received at z with frequency 
f in the range !-'o[(l + z)/{l + zs)] < v < i^o, where uq 
is the resonance line frequency, will have passed through 
the resonance line frequency en route to the absorber. To 
estimate the optical depth of the IGM at these frequencies, 
it is convenient to approximate the Voigt line profile as 



4>v{a,x) 



-1/2 / 2n 

TT exp(— a; ) 



l^l > X„ 



(Al) 



where Xm is the value of a; = (v — vo)/ Avd at which the two 
approximations match. For T = 10 K, 100 K and 1000 K, 
the values of Xm are 2.58, 2.83 and 3.05, respectively. 

A photon of frequency v a.t z will have x = iXm at the 
redshift zi^\ given by 



u 1 



,(±) 



1^0 1 + z 



= l±x, 



(A2) 



where b = (2fcBTIGM/T7^H)^''^ is the Doppler parameter 
of the IGM at temperature Tigm, assumed constant. (In 
principle, a turbulent velocity component may be added in 



quadrature.) The optical depth in Eq. ([T| may then be ex- 
pressed as a sum of three contributions, corresponding to 
the ranges x < —Xm, —Xm < x < Xm and x > Xm, 
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The first and third terms are similar to Eq. ([2]). The dom- 
inant term is the second, representing scattering through 
the line core. Setting u — [{^{1 + z')/(l + z) ~ uq) / Aud]'^ , 
and assuming b/c <^ 1 within the integrand, gives for this 
contribution 



1 aen(O) .3/2 1/2 3/2 .1 2 . 

8000h-\l + zf/^y-'/^'yi^,xl), 



(A4) 



where y = v/vo, 7(i, a;) = duu*~^ exp{—u) is an incom- 
plete gamma function, and Qm = 0.3 and Qth^ — 0.022 
are assumed. This will in general ensure that essentially 
no fiux from the source is received in the frequency range 
1^0(1 + •z)/(l + Zs) < V < vq. This is the Gunn-Peterson 
effect. 



APPENDIX B: SLAB SOLUTION 

The diffusion equation for the transfer of radiation through 
a static slab in t he Eddington approximation is given by 
iHarringtonl (|l973l ) as 



d nx{T,a) ^ a^n^(T, a) 



-3(j}v{a,x)^^ 
G{r) 



,1/2 



Snia) 



An 



(Bl) 



for a spectrally flat source Gij), in the limit cf>v{a,x) ~ 
a/{-Rx'^). Here, a is related to x through dx/da — 
{3/2)^^^(pv{a,,x), and r = T^/<l>v{a,x) is the mean opti- 
cal depth vertically through the slab. The photon number 
density n{x) is expressed as n^i^rj). The last line in Eq. (|B1|I 
follows from approximating 3(?!>v'(a, a;) as a Dirac 5- function, 
Q^^^Soicr) where the coefficient preserves the normalization. 
The boundary conditions assumed are 



dnx{r, a) 



dr 



= T^<Pv{a,x)n^{B,a) 



and 



ri^(r, a) for 



±00. 



(B2) 



(B3) 



For a uniform ly distributed source of unit strength, 
IHarringtonl ^97^) takes G = 1/(2B), where -B < t < B 
describes the vertical extent of the slab. In terms of the total 
optical depth at line centre, tq = 2_B(/>v(a, 0) = 2B/n^^^, the 
source becomes G = tt~^^^/to. The solution at the centre of 
the slab is 
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6l/2 



exp 



_2l/2 

V 3 / aro 



where the function S{z) = z - ^-^/S + 2''/5 - + . 

< 2 < 1. 

The diffusion approximation Eq. (|14p transforms to 



(B4) 
z for 



d'^n^icr) 1/2 4 r / N 

= -6 



(B5) 



for a source S{x) — ASo{x) and adopting D{a,x) = 
cj)v{a,x). Identifying A = G/(47r) = l/(47r^/Vo) repro- 
duces Eq. (IB1|I with the r-dependence suppressed. Express- 
ing Tcsc = Xf{aT(iY^^ , the solution Eq. (I18|l becomes 



127rV2 



(B6) 



for a;inj = 0. Noting 5(1) ~ 0.92, equating the values of 
Eqs. ([Bill and jBe]) at x = gives x, ~ 0.92. The two profiles 
are very similar, with Eq. (|B4|) forming an exponential tail 
with a cutoff at x > 0.87(aro)^/^ ~ 0.95a;osc. 
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